Biorthogonal Wavelet Based Algebraic Multigrid Preconditioners for Large Sparse Linear Systems
نویسندگان
چکیده
In this article algebraic multigrid as preconditioners are designed, with biorthogonal wavelets, as intergrid operators for the Krylov subspace iterative methods. Construction of hierarchy of matrices in algebraic multigrid context is based on lowpass filter version of Wavelet Transform. The robustness and efficiency of this new approach is tested by applying it to large sparse, unsymmetric and ill-conditioned matrices from Tim Davis collection of sparse matrices. Proposed preconditioners have potential in reducing cputime, operator complexity and storage space of algebraic multigrid V-cycle and meet the desired accuracy of solution compared with that of orthogonal wavelets.
منابع مشابه
A Novel Aggregation Method based on Graph Matching for Algebraic MultiGrid Preconditioning of Sparse Linear Systems
Multilevel techniques are very effective tools for preconditioning iterative Krylov methods in the solution of sparse linear systems; among them, Algebraic MultiGrid (AMG) are widely employed variants. In [2, 4] it is shown how parallel smoothed aggregation techniques can be used in combination with domain decomposition Schwarz preconditioners to obtain AMG preconditioners; the effectiveness of...
متن کاملNumerical study of the performance of preconditioners based on algebraic multigrid method and approximate sparse inverses
Application of algebraic multigrid method and approximate sparse inverses are applied as preconditioners for large algebraic systems arising in approximation of diffusion-reaction problems in 3-dimensional complex domains. Here we report the results of numerical experiments when using highly graded and locally refined meshes for problems with non-homogeneous and anisotropic coefficients that ha...
متن کاملForward Solving in Electrical Impedance Tomography with Algebraic Multigrid Wavelet Based Preconditioners
1. Abstract Electrical Impedance Tomography is a soft–field tomography modality, where image reconstruction is formulated as a non-linear least-squares model fitting problem. The NewtonRahson scheme is used for actually reconstructing the image, and this involves three main steps: forward solving, computation of the Jacobian, and the computation of the conductivity update. Forward solving relie...
متن کاملParallel PCG algorithms for voxel FEM elasticity systems
The presented comparative analysis concerns two parallel iterative solvers for large-scale linear systems related to μFEM simulation of human bones. The benchmark problems represent the strongly heterogeneous structure of real bone specimens. The voxel data are obtained by a high resolution computer tomography. Non-conforming RannacherTurek finite elements are used for discretization of the con...
متن کاملFunction-based Algebraic Multigrid method for the 3D Poisson problem on structured meshes
Multilevel methods, such as Geometric and Algebraic Multigrid, Algebraic Multilevel Iteration, Domain Decomposition-type methods have been shown to be the methods of choice for solving linear systems of equations, arising in many areas of Scientific Computing. The methods, in particular the multigrid methods, have been efficiently implemented in serial and parallel and are available via many sc...
متن کامل